Orientational orders of small anisotropic molecules confined in slit pores 
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We have studied phase behavior of hard gaussian overlap molecules with small anisotropic parameter 
confined in two plane parallel structureless hard walls. Our investigation based on standard constant- 
NPT Monte Carlo molecular simulation led us to some interesting findings. For small anisotropic 
molecules the nematic phase is instable in bulk, while, if the distance between the walls is small 
enough, an orientation-ordered phase can form. This result indicates that the required molecular 
elongation forming liquid-crystal phases is smaller in confinement than that in bulk. Considering the 
value of the elongation of molecules, the computed result inplies that small molecule liquid crystals 
may exist in confinement. 

PACS numbers: 61.30. Pq, 61.30.Cz, 61.20.Ja, 64.70.Md 
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I. INTRODUCTION 

There is a great number of molecular simulations in 
the literature on the properties of liquid crystals. Among 
these simulations, liquid crystal molecules are described 
using anisotropic models, such as the hard gaussian over- 
lap (HGO) models [1-3], the similar hard spherocylin- 
der [4-6] and the hard spheroidal models [7], and the 
wider used Gay-Berne (GB) models [8-11]. Although the 
detail results depend on the models used, some general 
conculsions have been found: for example it is generally 
confirmed that liquid-crystal phases can form only when 
the anisotropic parameters of these models are greater 
than some critical values. Otherwise (i.e. for molecules 
with small anisotropic parameters), as one increases the 
pressure or decreases the temperature, these systems will 
freeze or crystallize before forming liquid-crystal phases. 
Many recent works focus on the macroscopic properties 
and phase transitions of models with large anisotropy 
parameters. For example, based on molecular simulation 
and statistical theories, global phase diagrams of hard 
potential models [5,6] and GB models [8] in bulk have 
been developed. 

On the other hand, confinement in nanometer scale 
is found to induce phase transitions not observed in bulk 
systems and shifts in phase transitions [12-14]. Recently, 
Gelb et al. [12] systematically reviewed research trend on 
confined simple fluids based on statistical physics, molec- 
ular simulations, and carried out some experimental in- 
vestigations. The properties of confined liquid crystals 
also attracted wide interests on both sciences and tech- 
nologies, such as the surface anchoring effects of liquid 
crystals [15]. More recently, Gruhn and Schoen [10,16] 
studied the microscopic structure of confined liquid crys- 
tals based on GB models. So far, these researches on the 
bulk and confined properties of liquid crystals only fo- 
cus on systems composed of large anisotropic molecules 
which model the usual liquid-crystal molecules. To the 
best of our knowledge, properties of small anisotropic 
molecules have shown no interests, probably becuase they 



are believed not to form orientation-ordered phases in 
bulk [17]. 

But there is a need to address the following question: 
Are liquid crystal phases in confinement stable even in 
the small-anisotropic-molecule systems? According to 
the knowledge on confined simple fluids, although these 
systems cannot form orientation-ordered phases in bulk, 
confinement induce new phases when these molecules are 
confined in nanometer-sized pores. The answer to the 
above question is important not only for the sake of prop- 
erties of confined fluids, but implies the possibility of ex- 
isting small molecule liquid crystals. 

In this paper, we study phase behavior of molecules 
confined between two plane parallel structureless hard 
walls (i. e. slit pore), we investigate orientation-ordered 
phases and their transitions on a few kinds of molecules 
with different anisotropic parameters (molecular elonga- 
tion), and show that small anisotropic molecules can form 
stable orientation-ordered phase (nemetic phase) in con- 
finement, though in bulk the phase does not exist. Our 
results indicate that the required molecular elongation 
which provoke stable liquid-crystal phase in HGO sys- 
tems decreases with confining conditions. Therefore we 
suspect that small molecules might form liquid crystals 
in confinement. Some related discussions are presented 
in Sec. IV. Simulation results are shown in Sec. III. We 
begin, however, in Sec. II with a description of the HGO 
model and the presentation of our simulations. 



II. POTENTIAL MODEL AND 
COMPUTATIONAL DETAILS 

Due to a general belief that repulsive interactions are 
key roles to determine liquid cyrstal structures, some 
hard nonspherical models [1-7] are used to study proper- 
ties of liquid crystals. Among them, one of the most sim- 
plest models is the HGO model [1-3,17]. In the model, 
interaction between a pair of molecules i and j is given 
by 
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where, fjj = r^/r^ is a unit vector along the line join- 
ing the centers-of-mass of the molecules, in and iij are 
unit vectors along the principal axes of the molecules, a 
is dependent on the molecular orientation and explicitly 
given by 

ainj^Uj) = a Q {l -[ 1 + 
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Here, do is constant, which corresponds to the width of 
molecules, X is a measure of the nonsphericity of the 
molecule defined as X — (k 2 — l)/(k 2 + 1) with k being 
the length-to-width ratio of the molecule. 

We use the standard constant-pressure (const-NPT) 
Monte Carlo simulations to obtain the equation of state 
P{p) ( for hard potential, temperature T is not present 
explicitly ) . The simulations are performed with N = 108 
or 256 for the molecular elongations k being between 2.0 
and 3.0. For k > 3.0, a great number of simulations have 
shown that there is isotropic-nematic phase transition in 
the bulk [1-3], but for k < 2.0, the nematic phase is insta- 
ble in bulk up to a typical high density of liquids [17]. In 
simulations of bulk systems, the simulating box is cubic 
and equally fluctuates in three directions, cubic periodic 
boundary conditions are used. In confinement, the dis- 
tance between two walls (z direction) is fixed, the box 
equally fluctuates only in the two directions (x and y) 
and uses periodic boundary conditions. The simulations 
are organized in MC cycles, each MC cycle consisting (on 
average) of N trial translational and rotational molecu- 
lar displacements and one trial volume fluctuation. The 
maximum step length of each trial move is automatically 
chosen at each pressure for making the acceptable prob- 
ability fall between 0.4 and 0.6. The starting configu- 
ration is a fee lattice ( in confinement, the number of 
lattic in z direction may be smaller than that in other 
directions ), which melte at low pressure and equilibrate 
for 100,000 MC cycles. The system is slowly compressed 
in small pressure steps. For any given pressure, the sys- 
tem is typically equilibrated for 100,000 MC cycles and 
an average is taken over 100,000 additional MC cycles, 
then the final configuration is set as the starting config- 
uration of the next pressure. But 300,000 MC cycles are 
used near the isotropic-nematic transition. After arriv- 
ing at a high pressure point, in some cases, we slowly 
expand the system from this high pressure point to a low 
pressure point. In calculating average values (such as 
density, oricntational order parameter ), we sample par- 
tial equilibrated configurations. However, the obtained 
results (average values and their fluctuations) are almost 
the same as that using all equilibrated configurations. 

The orientational order parameter S is calculated in 
the simulations as the largest eigenvalue of the ordering 



Qa/3 tensor, defined in terms of the components of the 
unit vector Ui a along the principal axis of the molecules, 
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In the simulation of confined systems, a similar hard 
potential is set as the interaction between molecule and 
walls [10], thus 
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where z\ is the distance between center-of-mass of 
molecule i and wall k, and a w is given by 
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Here ij = (k 2 — l)/k 2 . We suppose the walls are at z = 
and z — D, respectively. Of course, the center-of-mass of 
the molecules must be between two walls. 

In the remainder of this paper, <r is set as unit of 
length, P/ksT is noted as P. 



III. RESULTS 

We first simulate the equation of state of HGO systems 
with k = 3.0 in bulk. For this case, isotropic-nematic 
(IN) phase transition had been found in previous simula- 
tions [1-3] . Based on free energy calculations of isotropic 
and nematic phases, the IN transition pressure was about 
5 for k = 3.0. In this paper, we do not try to deter- 
mine exactly IN transition point by using these free en- 
ergy methods, but only estimate the transition from their 
S — P curves. 

Our results are shown in Fig.(l), It is clear, the equa- 
tions of state exhibit a discontinuity at P ~ 5. In Fig. 
(1) (b), the behavior of the orientation order parameter 
S shows that a spontaneous IN phase transition occurs 
at the discontinuous point. The result in increasing and 
decreasing pressure (compressing and expanding simula- 
tions) are similar to each other and in agreement with 
previous simulations. 

To test the effects of confinement, we study the con- 
fined HGO systems with k = 3.0 in different D, the ob- 
tained results are shown in Fig. (2) and Fig. (3). This 
leads to our first interesting result: the IN transition 
point shift to lower pressure (i. e. density) direction as 
the distance between two walls D decreases. Usually, for 
small anisotropic molecules, the nematic phase is insta- 
ble, the reason can be understood. The system will freeze 
or crystallize before forming liquid-crystal phase while in- 
creasing the pressure or decreasing the temperature. So 
the shift at low pressure indicates that nematic phase 
form in confinement casilcr than in bulk. Therefore, we 
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suspect the required smallest k forming nematic phase 
is smaller in confinement than that in bulk. In Fig. (4) 
and Fig. (5), we show the equation of state of smaller 
anisotropic molecules (k = 2.0). Up to a typical high 
density of liquid ( packing fraction y is about 0.5), ne- 
matic phase is not found in bulk. Here y = pv m , where 
p is density of the liquid, v m = n/Qka^ is the volume of 
the molecule. The result is in agreement with previous 
simulations [17]. However, we do not find stable nematic 
phase in the system in confinement (the distance D be- 
tween walls is 3.5 and 4.0). Consider that k = 2.0 may 
be too small, we simulate another system with a little 
larger elongation (k = 2.2) and confined in a thiner pore 
(D = 2.5), the results are shown in Fig. (6) and Fig. (7). 
The system is compressed up to a typical high density ( 
y = 0.57 ), no phase transition is found in bulk. But for 
the confined system in slit pore with a part of D = 2.5, 
an obvious IN phase transition is found. In principle, to 
strictly determine the existence of liquid-crystal phase in 
bulk and in confinement, we should use free energy meth- 
ods to calculate the global phase diagram of the system 
including calculations of chemical potetials in isotropic 
phase, orientation order phase and crystallized phase, 
respectively. However, in this paper, we use a rougher 
method used by earlier Rigby [17]: compress the sys- 
tem to a typical high density (y is about 0.5), if there is 
no nematic phase in the density zone, the nematic phase 
does not exist, since in higher density zone, fluids usually 
freeze or crystalize. The result is a rough approximation, 
but we still can obtain a general qualitative conclusion: 
confinement will decrease the required molecular elonga- 
tion k to form a nematic phase. We believe that a ne- 
matic phase exists for a certain k in systems in suitable 
confinement, though the nematic phase of the systems 
is instable in bulk. Of course, it is important to obtain 
quantitative results on this effect, and we shall study 
global phase diagrams of more realistic small-anisotropic 
molecules in confinement. 



IV. CONCLUSION 

From the ongoing investigation, we can state that ori- 
entation ordered phase exists in bulk for some values of 
the molecular elongation k above a critical value k c ; for 
smaller value of k, the ordered phase becomes instable, 
but as these small k molecules are confined in very thin 
pores, the orientation ordered phase can be found. This 
result confirms our expectation that small-anisotropic 
molecules form liquid crystals in confinement. 

Usually, the elongation of liquid crystal molecules is 
about 4 — 5 and k of nonspheric inorganic molecules, such 
as CO2, CS2 etc., are only about 2. Hence, the usual 
small molecules cannot form liquid-crystal phases in bulk. 
But in confinement, according to our results, the required 
k forming liquid crystals will decrease, so it is interesting 
to study whether or not some small molecules will form 



liquid crystal phase in special confinement. Due to the 
wide interests in small molecule liquid crystals, it may be 
important to study phase behavior of real confined small 
molecules. On the other hand, complicated geometrical 
pores may induce more complex behavior in fluids [12]. 
Recently, Chiccoli et al. studied the properties of liq- 
uid crystals with dispersed polymer fibrils [18], where the 
polymer fibrils formed networks and served as complicate 
pores to confine liquid crystals. Considering the required 
sizes of pores to induce possible liquid crystal phases in 
small molecules may be very small (in molecule scale) , we 
shall study the effects of confinement in small anisotropic 
molecules with dispersed large anisotropic molecules. In 
the mixture, large anisotropic molecules may serve as 
complex pores to confine small molecules and induce 
possible orientation-ordered phases, which might induce 
characteristic of small-molecule liquid crystals. 
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Fig.1 X. Zhou et al. 108 HGOs with k=3 

FIG. 1. Isotropic-nematic phase transition in bulk is 
shown. Here, the number of molecules (N) in our simula- 
tions is 108. The anisotropic parameter k = 3. (a) Equation 
of state of the system P(p) is shown, (b) Orientation order 
parameter S versus the employed pressure. Shown simulation 
results include that an increasing pressure and a decreasing 
pressure process. 



Fig. 2 X. Zhou et al. 256 HGOs with k=3 confined in slit pores 
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FIG. 2. Equation of state P(p) of a HGO system in slit 
pores is shown. Here, N = 256, and k — 3. D is the width of 
the used slit pore in our Monte Carlo simulations. 



Fig. 3 X. Zhou et al. 256 HGOs with k=3 confined in slit pores 
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FIG. 3. Orientation order parameter S versus pressure P 
in slit pores with different D. A shift in IN transition is clearly 
shown. Here, iV = 256 and k = 3. 
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Fig.4, X. Zhou et al, N=256, and k=2.0 



Fig. 6, X. Zhou et al. N=108, k=2.2 
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FIG. 4. Equations of state P{p) of a HGO system in bulk 
and in confinement are shown. Here N = 256 and k — 2. The 
density of the system is increased to a typical high density of 
liquids from zero. 
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FIG. 6. Equations of state P(p) in a system in bulk and 
in confinement are shown, respectively. Here, N = 108 and 
k — 2.2. The densities of the system are increased to a typical 
high density from zero. In a thin slit pore (D — 2.5), the 
equation of state exhibits a discontinity. 



Fig. 5, X. Zhou et al. Confined HGOs with k=2.0, D=4.0 
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FIG. 5. Equations of state P(p) of systems are shown. Here 
N = 256 and 108. k = 2 and D = 4.0. Our result indicates no 
obvious sized effects in our simulations of confined systems. 



Fig. 7, X. Zhou et al, N=1 08, k=2.2 
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FIG. 7. Orientation order parameters S versus pressure P 
in a HGO system in bulk and in slit pores, where N = 108 
and k = 2.2. An IN transition is found in D — 2.5, but the 
transition does not exist in bulk. 



5 



